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X5 '. ABSTRACT 

D . 

Context. HD 61005 is a nearby young solar type star that shows a large infrared excess due to a debris disk. The disk has been recently 
, imaged from ground and space, with indications of several components. Some characteristics of the disk suggest the presence of 

planetary companions around the star, that remain undetected in deep adaptive optics imaging. 

Aims. For a better understanding of the system we aim to refine the determination of the stellar parameters, with emphasis on the 
stellar age and system orientation. 

, Methods. We used ASAS and Hipparcos photometry and FEROS spectra to determine the rotation period, radial and rotational 

velocity, chromospheric emission, effective temperature, and chemical composition. 

Results. We find no indication of any misalignment between the star rotation axis and the disk. The standard age calibrations applied 
' to several indicators yield an age close to that of the Pleiades (120 Myr); however the kinematic properties strongly support its 

CL(, membership in the younger (40 Myr) Argus association, which also includes the IC 2391 open cluster. Detailed comparison of 

the properties of HD 61005 and IC 2391 members shows that the characteristics of HD 61005 are compatible with membership to 
the Argus association, once its rather slow rotation is taken into account, because lithium and other age indicators are somewhat 
correlated with stellar rotation at a fixed age. We also identify systematic differences between the field and cluster population of the 
^ I Argus association, which are probably selection effects, so we suggest that additional members with slower rotation and lower activity 

^ level are waiting to be identified. 
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ff^ ^ 1 . Introduction the debris disks of the eariy type stars /3 Pic and FomaUiaut were 

, . ^ , , , , found to be accompanied by giant planets (Kalas et al. 120081 

. . HD 61005 IS a nearby (* 35 pc) solar type star thaUiosts a de- Lagrange et al. M and references therein). 
. bns disk discovered with Spitzer (Meyer et al. |2006). The in- 

O ■ frared excess with respect to the stellar photosphere at 24 fim , ^he interpretation of the direct imaging results in terms of 

is the largest that any star shows in the FEPS sample (Meyer et Pl^n^tary mass limits (Buenzli et al. 2010) and of possible fu- 

al. pOa ). The disk was resolved with HST (Hines et al. m ^""'^ planetary detections heavily relies on the derived stellar age. 

Maness et al. 2009) and VLT (BuenzU et al.OTI observations; [i;"^^ ' '^'"l*^ !ff of HD 61005 to be 90 ± 40 

the latter were obtained as part of the NaCo Large Program for ^y^^'^'^^, '2°^°.' Ro^catagliata et al. (2009) re- 

, Giant Planet Imaging (ESO program 1 84.C-0567). These stud- PO^'^ ^0 and 1 35 Myr, respectively. Considenng the special char- 

U ies revealed two components, an external one shaped by the in- acteristics of this system the uncertainty in the current age de- 

■ ■ ■ teractions with the interstellar medium and a ring seen nearly termination, and the availability of new observational data, we 

edge-on. The ring has a semimajor axis of 61.25 + 0.85 AU Performed a new study of the properties of the star, with the 

and an eccentricity of 0.045 ± 0.015, and its center is charac- 8°^! improvmg the stellar age determination. We also intend 

terized by an offset from the star of 2.75 ± 0.85 AU (BuenzU et '[^"J^ onentation of the star rotation axis and compare it 

al. 12010). Furthermore, Fitzgerald et al. (2010) report prelimi- ^^^8^-°" onentation of the debns ring (i = 84.3 ± 1 , 

nai-y evidence of a second, inner belt from 10 fim observations, ^ueji^li 2010) This is of special relevance considering the 

whose projection might be misaligned by about 20° with respect ^Jgh frequency of close-in planets whose normals to the orbital 

to that of the ring. These features might be explained by the pres- P^^"*^ ^o^ind to be misahgned with respectto the stellar rota- 

ence of planetary companion(s). Indeed similar characteristics in ^""^ al.EOEt. Recently, Lai et 

al. ( 20 1 ) have proposed that the interaction between the star and 

* Based on observations collected at La Silla and Paranal the disk might produce misahgned and even rettograde configu- 

Observatory, ESO (Chile). Programs 072.A-9006(A), 074.A-9020(A), rations. Systems with resolved disks and late spectral type that 

076.A-9013(A). Based on the All Sky Automated Survey (ASAS) pho- allow determination of the rotation period through photomet- 

tometric data. ric monitoring, such as HD 61005, are suitable targets for test- 



2 



S. Desidera et al.: The debris disk host star HD 61005: a member of the Argus Association? 



ing this hypothesis (Watson et al. 120101) . An alternative mech- 
anism that produces misaligned systems that might be at work 
for HD61005 is to add material from the interstellar medium 
onto a circumstellar disk, which alters the inclination of the disk 
(Moeckel & Throopl20T0l). 

Our paper is structured as follows. In Sect. |2]we describe the 
photometric data and the determination of the rotation period; in 
Sect. |3]we present the spectroscopic data and the data reduction 
procedures, and we describe the spectroscopic analysis we per- 
formed to derive radial and projected rotation velocity, chromo- 
spheric emission, effective temperature, and chemical composi- 
tion; in Sect. H) using the results of the previous sections as well 
as additional literature data, we discuss system properties such 
as orientation, kinematics, and binarity; in Sect.|5]we derive the 
stellar age and discuss the membership to the Argus association, 
and in Sect. |6] we summarize our conclusions. 

2. Rotation period 

We retrieved photometric time series of HD 61005 from both 
the ASAS All Sky Automatic Survey archive (Pojmanski 2002 1 
(data collected from November 2000 to November 2009) and the 
Hipparcos and Tycho epoch photometry archives (data collected 
from November 1989 to February 1993). 

For Hipparcos and Tycho photometry, we selected 240 and 
183 measurements, respectively, quoted as quality A or B. After 
removing very few outliers, we were left with 236 and 180 mea- 
surements, respectively. From Hipparcos data we derived a mean 
magnitude of V=8.22 and a standard deviation nns-0.024 mag, 
whereas the data precision was cr=0.015 mag. From Tycho data 
we derived a mean magnitude V=8.21 and a standard deviation 
rms=0.13 mag. However, the data precision was a much poorer 
(T=0. 1 1 mag. The Hipparcos and Tycho magnitudes were trans- 
formed into the Johnson standard system according to the equa- 
tions in Ferryman et al. (1 19971 1. 

For ASAS photometry, we retrieved 617 useful measure- 
ments (quality A-C). Among five different apertures, we found 
that a 4-pixel (~30" radius) aperture photometry extraction pro- 
vided the highest photometric accuracy. After removing outliers 
(by applying a 3cr threshold) we were left with 588 magnitude 
values useful for subsequent period search. The mean magni- 
tude was V=8.22 mag with a dispersion of rms=0.028 mag and 
an average photometric precision of cr=0.03 mag. The typical 
time sampling of ASAS observations was one observation every 
two to three days. 

2.1. Search for rotation period 

A visual inspection of the HD 61005 photometric time series 
clearly shows this star is variable. Based on its late spectral 
type, we expect this variability to arise from brightness inho- 
mogeneities (starspots) that are carried in and out of view by 
the star's rotation, unevenly distributed across the stellar photo- 
sphere. 

To search for significant periodicities in the HD 6 1 005 ASAS 
time series coming from stellar rotation, we sectioned the com- 
plete series into 14 time intervals never exceeding about two to 
three months. This approach has the twofold advantage of mini- 
mizing possible phase shifts introduced by the growth and decay 
of active regions, and of allowing us to reveal possible temporal 
variations in the rotation period that may be related to surface 
differential rotation. 

The Lomb-Scargle periodogram (Scargle 1982] was com- 
puted on each segment, as well on the complete time series. In 



the latter case, provided that the light curve remains quite stable 
(as proved to be the case of HD 61005), we can detect the rota- 
tion period with a much higher precision and confidence level. 
A more detailed description of the data preparation and period 
search (determination of false alarm probability and period un- 
certainty) can be found in Messina et al. (I2010al l. 

2.2. Results 

We found that HD 61005 has a rotation period of f = 5.04+0.04 
days. The same rotation period is found in the ASAS, Hipparcos, 
and Tycho photometry. In Table [1] we Ust the results of our 
period search. We list the rotation period, its uncertainty, the 
peak-to-peak light curve amplitude, the photometric precision, 
the JD time interval in which the period search was performed, 
the number of measurements, the normalized peak power in the 
Lomb-Scargle periodogram (at 99% confidence level, as com- 
puted from simulation), and the normalized power derived from 
real data, the brightest magnitude, the source of photometry, and 
the results derived from the time series without sectioning. 

All the rotation period determinations are consistent with 
each other within the computed errors. The photometric period 
variations allow derivation of a lower limit to the amplitude of 
the surface differential rotation (see Messina & Guinan 2003). In 
the case of HD 61005, these variations are less than 1%, whereas 
the average period uncertainty is AP = 3%. Therefore, based on 
the available photometry we can exclude the presence of signif- 
icant (> 1% at 3cr level) surface differential rotation. 

In Fig. [1] we plot the complete ASAS magnitude time series 
(top panel) together with a selection of time segments. In the left 
panels, we plot the timeseries magnitudes along with an over- 
plotted sinusoid (solid curve) with a P=5.04d rotation period. In 
the middle panels, we plot the periodograms with indication of 
the 99% confidence level (horizontal dashed line) and mark the 
power peak cotTesponding to the rotation period. Finally, in the 
right panels we plot the phased light curves using the P=5.04d 
period and the first JD of the time segment as the initial epoch. 
Here again the solid curve is the fitting sinusoid. In the third 
panel from top, we notice that there may be time segments where 
the light curve has a clear rotational modulation, although the 
power peak does not reach the 99% confidence level. 

2.3. Long-term piiotometric variability 

The Lomb-Scargle periodogram analysis performed on the 
joined ASAS and Hipparcos datasets also reveals a periodic 
long-term variability with a period of P(;ycie=3. 4+0.1 yr and 
a confidence level larger than 99%. The amplitude appears to 
change from cycle to cycle, as generally observed in cyclic spot- 
ted stars, and reaches a maximum value of about A^yde -0.08 
mag. In Fig.|2]we plot the Hipparcos + ASAS V band magni- 
tude time series with an overplotted fitted sinusoid of P=3.4yr 
cycle. Similar results are obtained when the periodogram anal- 
ysis is carried out on only the ASAS dataset. The Tycho mag- 
nitudes are excluded from this analysis, since their uncertainties 
are much greater than contemporary Hipparcos magnitudes. 

The long-term variation probably represents a starspot cy- 
cle consisting of a quasi-periodic change of the star's spot- 
tedness similar to the 11 -yr sunspot cycle. It is very interest- 
ing to note that rotation period, cycle period, and cycle ampli- 
tude of HD61005 are all very similar to the values observed 
in the Pleiades-age spotted star DX Leo that has Prot=5.42d, 
Pcycie=3.1±0.05, and Acycie -0.12 mag (Messina & Guinan 
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Fig. 1. The complete ASAS magnitude time series (top panel), 
together with a selection of time segments. In the left panels, 
we plot the time series magnitudes with an overplotted sinusoid 
(solid curve) of P=5.04d rotation period. In the middle panels, 
we plot the periodograms with indication of the 99% confidence 
level (horizontal dashed line) and the power peak coiTesponding 
to the rotation period. Finally, in the right panels we plot the 
phased light curves using the P=5.04d period and the first JD of 
the time segment as initial epoch. Here again the solid curve is 
the fitting sinusoid. 



I2002I I. Both stars are well placed on the active branch observed 
by Saar & Brandenburg (11999) for the distribution of the activ- 
ity cycle frequency (normalized to the rotation frequency) versus 
inverse Rossby number for active stars. 

Unfortunately, the B-V color from Tycho, which has quite 
a large average uncertainty of crB_v=0.15 mag, is the only 
one available to investigate the expected color variations of 
HD 61005 originating from magnetic activity. We averaged the 
Tycho series of B-V measurements (after removing points with 
uncertainties above 0.20 mag) using time bins of 200 days in or- 
der to explore possible color variations. However, the inferred 
long-term B-V variations are found to be a factor 2.5 smaller 
than the average uncertainty. A Pearson linear correlation test 
shows such long-term B-V variations to be uncorrelated at a sig- 
nificante level a > 50% (cfr. Bevington ll969l l. In the bottom 



Fig. 2. Top panels: the periodogram (left) and the window func- 
tion (right) obtained from the Lomb-Scargle analysis of the 
ASAS plus Hipparcos time series. Middle panel: Hipparcos + 
ASAS V band magnitude timeseries with an overplotted fitted 
sinusoid of P=3.4yr cycle. Bottom panel: V versus B-V mea- 
surements from Tycho (individual measurements as well as the 
binned values). 



panel of Fig. |2l we plot the V versus B-V measurements from 
Tycho, along with the binned values. We notice that both individ- 
ual measurements and average values show the tendency of the 
star to be bluer when fainter. Owing to the very poor quality of 
the available photometry, this evidence must be taken with great 
caution. If the observed trend is real, it may suggest that besides 
cool starspots, the observed photometric variability of HD 6 1 005 
also arises from bright faculae (see, e.g., Messina ,2008; . 

2.4. Photometric color 

Beside photometric colors in Hipparcos catalogs (B-V=0.734, 
V-I=0.78, from Tycho and Star Mapper respectively), high- 
quality broad band photometry was performed by Menzies et 
al. (1990). It does not fully sample stellar variability (mean 
of 8 measurements); however, the V magnitude (8.225+0.008) 
agrees very well with the long-term ASAS and Hipparcos 
magnitudes (see above). Therefore, in the following we adopt 
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the B-V=0.751 ±0.003 and V-I=0.805+0.005 from Menzies et 
al. (1 19901 1. Near-infrared magnitudes are available from 2MASS 
(J=6.905+0.026; H=6.578+0.046 Ks =6.458+0.024). We obtain 
V-Ks ^ 1.762 + 0.025. 



3. Determination of spectroscopic parameters 

We exploited 12 spectra of HD 61005 taken in five nights be- 
tween December 2003 and March 2006 using FEROS. These 
data were part of the RV survey for planets around young stars 
by Setiawan et al. 62008^ . The FEROS spectra extend between 
3600 A and 9200 A with a resolving power R = 48000 (Kaufer 
et al. I1999l l. Three spectra were acquired using the object-sky 
set-up and the remaining using object-calibration mode (see 
Table |2ll. The data were reduced using a modified version of the 
FEROS-DRS pipeline (running under the ESO-MIDAS context 
FEROS), which yields the wavelength-calibrated, merged, nor- 
malized spectrum. The reduction steps were the following: bias 
subtraction and bad-column masking; definition of the echelle 
orders on flat-field frames; subtraction of the background diffuse 
light; order extraction; order by order flat fielding; determination 
of wavelength-dispersion solution by means of ThAr calibration- 
lamp exposures; order-by-order normalization, rebinning to a 
linear wavelength-scale with barycentric correction; and merg- 
ing of the echelle orders. 

3.1. Projected rotational velocity and radial velocity 

The radial-velocity (RV) and vsin / were determined by cross- 
correlating the object spectra with a low-vsin/ template spec- 
trum (i.e. a solar spectrum) obtained with the same instrument 
and reduced in the same way. The normalized target spectrum 
was preliminarily rebinned to a logarithmic wavelength scale 
(A In /I = Av/c = 5 X 10"^) and split into six wavelength- 
ranges that are free of emission lines and of telluric absorptions 
(cf. Esposito et al. I2007I I. A Gaussian was fitted to the peak of 
the cross-correlation function (CCF), computed in each of the 
six distinct spectral ranges. The resulting RV values were av- 
eraged with weights (Schisano et al. 2009). The projected rota- 
tional velocity v sin / was derived from the FWHM of the cross- 
correlation peak in each of the six ranges above. The relation 
FWHM-vsini was obtained by convolving the template spectrum 
with rotational profiles of vsini from 1 to 100 km/s. Our mea- 
surements are listed in Table |2] 

In Table [3] we collect the measurements of v sin / available 
for HD 61005. All the determinations except the one by White et 
al. (I2007I I. which is based on data with lower spectral resolution, 
show fairly good agreement. In the following analysis we adopt 
our own determination (v sin / = 8.2 ± 0.5 km/s). 

3.2. Teff determination 

The determination of the effective temperature was obtained 
through the method of equivalent widths (EW) ratios of spectral 
absorption lines, using the calibration for FGK dwarf stars by 
Sousa et al. (2010). The EWs were measured using the ARESQ 
automatic code (Sousa et al. 12007b . In order to test the perfor- 
mances of ARES on spectra with different rotational broaden- 
ing, we used a FEROS spectrum of the Sun broadened artifi- 
cially with rotation profiles of increasing v sin /. The resulting 
values of effective temperature remain consistent (within 50K) 
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Fig. 3. Effective temperature derived with ARES on Sun spectra 
broadened artificially to various projected rotational velocities. 
Results are robust for v sin / up to 18 km/s. 

with the accepted value for the Sun up to a rotational velocity of 
about 18 km/s (Fig. |3]l. For higher rotation velocity, the number 
of measured line-ratios drops drastically with correspondingly 
larger errors due to increased line blending. The temperature de- 
terminations become unreliable for v sin / approaching 30 km/s. 

The value of Teff determined by ARES (5550 K) is con- 
sistent with those derived from excitation equilibrium as part 
of the abundance analysis (5500+60 K; Sect. [33] ). from B-V, 
V-I and V-K colors using the calibration by Casagrande et 
al.|20l0!(5480 K), the relationship between W(NalD) vs. Teff 
by Tripicchio et al. (1997) (5550 K), and the previous determi- 
nation by Hines et al. (2007) (5456 K from SED fitting using 
Kurucz models). 

3.3. Abundance analysis 

There was no previous abundance determination based on high 
resolution spectra. Gray et al. (2006 1 derived [M/H]= -0.06 
from modeling of their low resolution spectrum. The lower value 
derived by Nordstrom et al. (.2004) from Stromgren photometry 
([Fe/H]= -0.19) is biased by the effect of stellar activity (see 
Favata et al. [T997] l. 

We performed elemental abundance analysis for HD 61005, 
measuring abundances for several iron-peak and (T-elements, 
namely Fe, Ni, Na, Si, Ca, Ti, and the .s-process element Ba. 
The chemical analysis was carried out following the same ap- 
proach and methodology as in D'Orazi & Randich (2009, here- 
after DR09). The Ba content was instead derived as in D'Orazi 
et al. (I2009al) . taking isotopic splitting and hyperfine structure 
into account for both Ba ii features at 5853 A and 6496 A (see 
D'Orazi et al. 2009a for further details on this issue). 

Our procedure can be briefly summarized as follows. We em- 
ployed the line list provided in DR09, deriving LTE force-fitting 
abundances with the driver abfind in MOOG (Sneden 19731 2002 
version) and the Kurucz (119931 ) set of model atmospheres. The 
EWs were measured with the IRAF task SPLOT, performing 
a Gaussian fit to line profile. All the features with EW values 
higher than ~150 mA were discarded to avoid spectral lines 
close to the saturation regime of the curve of growth. 

Stellar atmospheric parameters were spectroscopically opti- 
mized: effective temperature (Teff) and microturbulence (^) val- 
ues were obtained zeroing the slope of log n(Fe i) with the 
excitation potential (x) and with line strengths, respectively. 
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Concerning gravity, we imposed the ionization equilibrium con- 
dition, i.e., log n(Fei)-log n(Feii)<0.05 dex, to derive the log g 
value. 

We obtained Teff=5500+60K, ^=1.00+0.25 km/s, and 
log g=4.5+0.2 as final adopted parameters. We refer the reader 
to DR09 for a detailed discussion on error estimates; here we just 
mention that internal errors, due to EW measurements and atmo- 
spheric parameters, are not larger than 0.08 dex. As to systematic 
(external) errors, our analysis (with the same method, code, line 
list) of the Hyades star VB 187 indicates that no major uncer- 
tainties affect our abundance estimates (see DR09). In Table |4] 
[Fe/H] and [X/Fe] ratios are reported for HD 61005 along with 
internal uncertainties due to EW values. The [X/H] ratios were 
computed by assuming the ones derived by DR09 as solar abun- 
dances. As a comparison, average abundances from seven stars 
in the open cluster IC 2391 are given and the standard devia- 
tion from the mean (rms) is provided in coiTespondence of each 
[X/Fe] mean ratio. 



3.4. Lithium 

The EW of the Li 6708 A resonance line was determined by 
Gaussian fitting with the ARES code, (see Table |2|. Our result 
(EW=171 mA) agrees closely with the determinations available 
in the literature (Table |3}- 

The Li i feature at 6707.78 A is an unresolved doublet so 
we conducted a spectral synthesis analysis to derive a reliable Li 
abundance. Using the driver synth in MOOG and Kurucz model 
atmospheres (as previously done for our abundance analysis, see 
Sect. 13.3b . we computed synthetic spectra in a wavelength win- 
dow of ~ 15 A, from 6700 A to 6715 A, and employed the 
line list given in D'Orazi et al. (I2009bl l. As a first step, we opti- 
mized the line list by changing the oscillator strengths, i.e. log gf 
values, when necessary, to obtain the best agreement between 
our FEROS solar spectrum and the standard iron abundance by 
Anders & Grevesse (.1989) . To perform the comparison between 
the observed spectrum of HD 61005 and the synthetic ones, we 
convolved the latter with both a Gaussian profile at our reso- 
lution of R~48000 (FWHM=0.14 A) and a rotational profile, 
taking the limb-darkening coefficient into account. 

Concerning atmospheric parameters, for consistency we 
adopted the ones derived from our EW abundance analysis, 
namely TefF=5500K, log g=4.5, ^=1.0 kms ', and [Fe/H]=0 
dex. As a by-product, the spectral synthesis allowed us to infer 
the projected rotational velocity and we obtained vsini=9.0+0.5 
kms"', in perfect agreement with other available estimates (see 
TableEJ. 

As shown in Fig. |4] all the observed features are reproduced 
very well by the synthetic profiles, which simultaneously fit the 
Fe I fines (6703, 6704, and 67 10 A) and the Li i doublet. The best 
agreement between observed and synthetic spectra is achieved 
for a Li abundance of log n(Li)=2.85 (righthand panel of Fig.|4]i. 

Uncertainties related to the best-fit determination and to the 
adopted Teff are 0.05 dex and 0.07 dex, respectively (note that Li 
abundance is only scarcely affected by other stellar parameters, 
i.e., log g and ^); to estimate the total error, we quadratically 
summed both these contributions, given their independence, and 
we found 0.09 dex. In passing, we also recall that uncertain- 
ties due to non LTE (NLTE) and/or three-dimensional model at- 
mospheres can be neglected, since they cancel each other (see 
Asplund & LindlMnifor details). 




6704 6706 6708 6710 671S 6706 6707 6708 6709 

Wavelength (A) Wavelength (A) 



Fig. 4. Left panel: Spectral synthesis of the region around the 
Li 6708 A doublet. Right panel: details of the Li doublet. 
Overplotted is our best-fit Li abundance (log n(Li)=2.85; red 
short-dashed curve) while the dotted and long-dashed curves 
represent log n(Li)=2.75 and log n(Li)=2.95, respectively. 



3.5. Chromospineric emission 

Chromospheric emission in Call H&K lines of HD 61005 was 
determined by measuring the S index on FEROS spectra as in 
Desidera et al. (12006b . The S index was also measured on 73 
FEROS spectra of 41 stars from the lists of Baliunas et al. (11995b 
and Wright et al. (12004b . reduced in the same way as those of 
HD 61005, to calibrate our instrumental S index onto the stan- 
dard M. Wilson scale. The dispersion of the residuals of the cal- 
ibration is 0.026, most likely dominated by the intrinsic vari- 
ability of chromospheric emission. The variability of the S in- 
dex measured on 15 spectra of the chromospheiically quiet star 
T Ceti is 0.0034 (2.0%). Full details of the calibration will be 
provided in a forthcoming paper. There are other four measure- 
ments of H&K emission of HD 61005 in the literature listed in 
Table [3] The average of our nightly averages and the literature 
determinations (excluding those by Gray et al. 2006 and White 
et al. 2007 because of lower accuracy) yields a mean S index of 
0.502 (rms 0.017), which corresponds to log - -4.310. The 
expected rotational period for a star with such activity level and 
color using the calibration by Mamajek & Hillenbrand (12008b is 
5. Old, very close to the observed one (5.04d). 



4. System parameters 

4. 1 . Photometric and trigonometric distances 

The trigonometric parallax in Van Leuveen (2007^1 is 28 .29+0.85 
mas. The position in the color-magnitude diagram agrees with 
those of the members of other groups in the age range 30- 
100 Myr, while ages under 20 Myr (see Fig. |5]l are excluded. 
Estimated photometric distances (from B-V, V-I and V-K col- 
ors) using the sequences of AB Dor and Tucana associations are 
37 and 39 pc, respectively, a not significantly different from the 
trigonometric one (35.3+1.1 pc. Van Leuveen l2007b . which we 
adopt in the following. 
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Fig. 5. My vs B-V, V-I and V-K color magnitude diagrams for 
HD 61005 (red star), IC 2391 ope n cluste r (black filled circles; 
data summarized in Messina et al. I2010bl l. /3 Pic moving group 
(orange filled triangles, 10 Myr), Tucana, Carina and Columba 
associations (green filled squares, 30 Myr) and AB Dor moving 
group (blue filled, upside-down triangles, 70 Myr). For these as- 
sociations, only the stars with trigonometric parallaxes among 
the members selected by ToiTes et al. (2008 1 were considered. 
The position of HD 61005 is compatible with a broad range of 
ages, older than about 20 Myr. 

4.2. Teff, luminosity, and radius 

Using the absolute visual magnitude My = 5.478 and the bolo- 
metric coiTection -0.162 from Flower (1 19961 1, we derive a lumi- 
nosity L - 0.583 Lq. Coupled with the adopted effective tem- 
perature of Teff =5500 K (Sect. [32]), this gives a stellar radius of 
R = 0.840 + 0.038 Rq (compatible within errors with the deter- 
mination by Watson et al. l2010i R = 0.829 + 0.048). 

4.3. System inclination 

The measured rotational period of 5.04d coupled with the stel- 
lar radius of 0.840 + 0.038 Rq yields a rotational velocity of 
8.43 + 0.38 km/s. Coupled with our adopted projected rotational 
velocity of vsini = 8.2 + 0.5 km/s, this gives sin / = 0.97 + 0.09 
and i = 77°^ J^. The ring inclination is 84.3° + 1 (Buenzli et 
al. 1201 Oi l, with a corresponding sin / of 0.995. 

We conclude that both the star and the disk are viewed close 
to edge-on and that there are no indications of any misalign- 
ment. A similar conclusion was reached by Watson et al. (2010). 
However, their analysis is based on the expected rotation period 
from H&K emission and the Noyes et al. ( 1984) relation rather 
than the measured value, which was derived for the first time in 
our study. 



4.4. X-ray emission 

The star is present neither in the ROSAT all-sky survey (RASS) 
catalog nor in other X-ray catalogs. We note that there are no 
sources in the official RASS catalogs within 1.2 deg from HD 
61005. This is unexpected considering the young age inferred 
from other indicators and the distance from the Sun. However, 
in the paper by Wichmann et al. (I2003I I. an X-ray luminosity 
\ogLx - 29.2 for the star is quoted, which was based on a 
preliminary RASS catalog, which is older than the official one. 
The source is not in the official catalog most likely because it 
may have been discarded as spurious in the final RASS data 
screening (Wichmann 2010, private communication). Therefore, 
the original RASS data have been kindly revisited by Jiirgen 
Schmitt. The source is clearly recognizable by eye. The EXSASQ 
source detection algorithm gives a broad-band count rate of 
0. 176+0. 034cts/sec (32.21+6.27 counts in 272.4 sec exposure 
time, and vignetting correction factor of 1.493) and a hard- 
ness ratio practically equal to zero. Adopting the calibration by 
Hunsch et al. (1999) and the trigonometric distance, an X-ray 
luminosity log Lx = 29.3 ± 0.2 is derived. 

4.5. Radial velocity, binarity and planetary companions 

There are no indications of any RV variations with amplitudes 
larger than 1 km/s from the various literature data or from our 
own RV measurements based on FEROS data. We therefore ex- 
clude the star being is a spectroscopic binary. The RV curve by 
Setiawan et al. (I2008I I shows a significant scatter with peak-to- 
peak amplitude of about 150 m/s and a correlation with line pro- 
file. Therefore, the RV variability appears dominated by the line 
profile alterations caused by stellar activity, and the limits on 
close-in planets are rather coarse, excluding only massive plan- 
ets in close orbits. The deep imaging also rules out either stel- 
lar, brown dwarf, or even massive planetary companions at a 
projected separation larger than a few tens of AU (Buenzli et 
al. l2010l l. The similarity between short-term Hipparcos and his- 
torical proper motions also point to a single star. Overall, there 
are no indications of any massive companions at either short or 
wide separations. The eccentricity of the ring suggests a giant 
planet within a few tens of AU (Buenzli et al. 1201 Oi l. 

4.6. Kinematic parameters 

The space velocities of HD 61005 are U = -22.2 ± 0.6 km/s, 
V = -14.3 ± 0.3 km/s, and W = -4.1 + 0.2 km/s. These are very 
similar to those of the Argus association (-22.0 ± 0.3 - 14.4 ±1.3 
-5.0 + 1.3) as given in Torres et al. (12008) . The space position 
is X = -13.9 ± 0.4 PC Y = -32.3 ± 1.0 pc, andZ = -3.3 + 0.1 
pc, which places the star at the outskirts of the known members 
of the Argus association (Fig|6]l. 

The existence of the Argus association was first proposed by 
Torres et al. (2003) from SACY data, thanks to its peculiar U 
velocity. The young open cluster IC 2391 was shown to share 
the motion and the age of the proposed association. The occur- 
rence of groups of young field stars associated with IC 2391 
had been formerly proposed by Eggen (19911 and Makarov & 
Urban (2000); however, only a minor fraction of their proposed 
members is part of the Argus association according to Torres et 
al. (2008) . 



- EXSAS stands for Extended Scientific Analysis System to evaluate 
RASS data; see Zimmermann et al. ( 11994) for details 
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Table 5. Summary of age determination for HD 61005 
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Fig. 6. Upper panels: U,V, and W space velocities forHD 61005 
(red asterisk) and Argus members from the list of Torres et 
al. l2008l Square symbols refer to objects with trigonometric par- 
allaxes. Lower panels: space positions X, Y, Z. Symbols as in the 
upper panels. 



Torres et al. (|2008t give an updated list of members of this 
association separated as field stars and IC 2391 members (29 and 
35 stars, respectively, see also the review by Pettersson I2008I I. 
The estimated age of the association is 40 Myr, in the middle of 
the proposed ages for IC 2391 (30+5 Myr Staufferetal. 199Tus- 
ing isochrone fitting, 50+5 Myr Barrado y Navascues et al. .20041 
using the Li depletion boundary, 30-50 Myr Platais et al. l2007l l. 

The mean distance of the members in Torres et al. (20081 is 
106 pc, much farther than the 35 pc of HD 61005, but a few in- 
dividual members have comparable distance. The trigonometric 
distance to IC 2391 is 144.9 ± 2.5 pc (Van Leeuwen 2009), and 
it is expected that several field members at distances compara- 
ble to or larger than IC 2391 wait to be identified (assuming a 
roughly spherical spatial distribution). 

HD61005 is not included in the member list by Torres et 
al. (2008). This can be easily understood, because the SACY 
sample (Torres et al. 2006) is based on a cross-identification of 
the ROSAT All Sky Bright Source Catalog (Voges et al. IT999] ), 
of which HD 61005 is not included (see Sect.l4!tl). 



The young age clearly resulting from all the indicators de- 
scribed above and the very similar space velocity call for a de- 
tailed evaluation of the possible membership of HD 61005 to the 
Argus association. This is the goal of the next section. 



5. The age of HD61005 and the membership to 
Argus association 

In the previous section we measured or collected from the liter- 
ature the determination of several parameters that are known to 
depend on stellar age. We discuss here whether they are compat- 
ible with the Argus membership suggested by kinematic prop- 
erties or if they instead indicate a diff'erent age. Application of 
existing age calibrations is summarized in Table |5] The calibra- 
tion of lithium EW into age was built using literature measure- 
ments of nearby young associations and Pleiades and Hyades 
open clusters and will be presented elsewhere. 



Indicator Age Calibration 



Li 


113 


see text 


H&K 


148 


Mamajek & Hillenbrand ( 2008J 


H&K 


72 


Donahue (1993 ); Henry et al 1 19961) 


Xray 


108 


Mamajek & Hillenbrand (20083 


P 


186 


Mamajek & Hillenbrand {,200% 


P 


125 


Barnes (2007 J 



The ages in Table|5]are in general older than the accepted age 
of the Argus associatiorQ. At first sight, this argues against mem- 
bership. However it is known that most (possibly all) of the HD 
61005 age indicators have some correlation with stellar rotation, 
and show a wide dispersion in coeval (especially young) stel- 
lar groups and clusters. A one-to-one correspondence between 
color and rotation period is definitively reached by an age of 
about 500-600 Myr (see, e.g.. Collier Cameron et al. 2009 ). The 
dispersion in rotation rate for coeval systems might be linked 
to the disk lifetime (disk-locking scenario). The active debris 
disk of HD 61005 might be a relic of a long-lived primordial 
disk. This might have prevented the system reaching fast rotation 
rates, increasing at the same time the lithium depletion (Bouvier 
,2008 J . A more detailed discussion is therefore needed to see if 
there are any known Argus members that are slow rotators and 
present similar discrepancies in the age indicators to HD 61005. 
Furthermore, we also consider here additional membership cri- 
teria such as chemical composition. 



5.1. Rotation period 

The rotation period of HD 61005 was determined in Sect. |2] 
while those of Argus and IC2391 members are taken from 
Messina et al. (1201 Obi ). They classified the detection of rota- 
tion period as either confirmed, likely, or uncertain depending 
on the number of significant detections of the periodicity among 
the segments by which the photometric time series were di- 
vided. We use this classification in the following discussion. The 
four stai-s (CD-582194 and CD-621 197 in Argus, PMM351 and 
PMM 2182 in IC2391) for which the rotation period is inconsis- 
tent with the spectroscopic v sin / are not used in the analysis. 

The rotation period of HD 61005 compared to that of mem- 
bers of Argus (field) and IC2391 is shown in Fig. |7] The dis- 
tribution of rotation periods of IC 2391 stars appears different 
from that of Argus field stars, with more slow rotators in the 
cluster. This might be from the different selection criteria of the 
IC 2391 cluster and the Argus association field members. Argus 
field members are identified in the sample of the SACY sur- 
vey (ToiTes et al. 2006 ) and, by definition, have X-ray counter- 
parts in the ROSAT Bright Source Catalog (Voges et al. ,19991 . 
Conversely, the members of IC 2391 are selected from a variety 
of sources, including deep X-ray imaging of the central regions 
of the cluster (Patten & Simon 1996; Simon & Patten 19981 
Marino et al. 2005 ) or a combination of spectroscopic, astromet- 
ric and photometric criteria, without X-ray preselection (Platais 
et al. 2007). Indeed, several members identified by Platais et 
al. (2007) only have X-ray counterparts in the ROSAT Faint 
Source Catalog (Voges et al. 120001 ). and a few targets have no 



^ Weise et al. ( [20T0l estimated an age of 30 Myr for HD 61005 from 
Li EW. This is partially due to the younger ages assumed for the clusters 
adopted as calibrators. 
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Fig. 7. Upper panel: distribution of rotation periods of Argus as- 
sociation (blue squares) and IC 2391 (black circles) members 
versus B-V color. Confirmed and likely periods are plotted with 
filled symbols, uncertain periods with open symbols. Data are 
taken from Messina et al. (I2010bl l. HD 61005 is represented 
by the (red) star symbol. The dashed and dotted lines represent 
the upper bounds of the period distribution of Pleiades members 
(from Bames [2003l l. 



X-ray counterparts at all (see App. |A}. Since fast rotators have 
brighter coronal luminosities and the mean distance of Argus as- 
sociation members is 106 pc, we expect that the census of associ- 
ation members is biased toward the most active and fast-rotating 
stars. We therefore consider the IC 2391 cluster members as a 
better comparison sample than Argus field members to evaluate 
the possible link of HD 61005 with this association. 

The position of HD 61005 in the color-period plot of Fig. Q 
is on the margin of the IC 2391 period distribution, close to a few 
members (PMM665, PMM1373, PMM4362, and PMM6974) 
whose rotation period still need confirmation. HD 61005 appears 
to be along the so-called "interface sequence" that represents 
the upper bound of the period distribution. Here this sequence 
is represented by the gyro-isochrone from Barnes (t2003t and 
computed for the Pleiades age of 120 Myr. 

We further note that the photometric amplitude is within the 
distribution of Argus and IC 2391 members with similar rota- 
tional period (see right panel of Fig. 6 in Messina et al. l20I0bl l. 

5.2. Lithium 

The comparison of Li EW of HD 61005 with that of Argus as- 
sociation members (Torres et al. 12006) shows that HD 61005 
has an EW about 50 mA lower than that of members of similar 
color (Fig.[8]l. When one instead considers members of the open 
cluster IC 2391 (from Platais et al. 12007. Randich et al. ISUOTl 
Stauffer et al. 1997 ), HD 61005 falls on the edges of the distri- 
bution. The Li EW also fits within the distribution of the Pleiades 
observations (age about 120 Myr). Application of an age-EW Li- 



■o o, ' o 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 

B-V 



Fig. 8. Comparison of Li EW of HD 61005 (red filled star) 
with that of Argus association field stars (blue squares) and IC 
2391 open cluster stars (black circles) Filled symbols refer to 
stars with confirmed or likely rotational period in Messina et 
al. (2010b), and open symbols to stars with uncertain or unde- 
termined rotational period in Messina et al. (l2010bV The dashed 
curve is the cubic fit to the locus of IC 2391 members in the 
range B - y = 0.50 - 1.10 



color relation based on data of nearby young associations and 
Pleiades and Hyades open clusters yields 113 Myr. 

It appears that, in spite of the association between IC 2391 
and the Argus field stars, IC 2391 members have on average a 
slightly smaller lithium content. At first sight this can be consid- 
ered a signature of a small age difference between cluster and 
field members. However, considering the possible dependence 
of lithium abundance on rotational period or rotational velocity 
found by Soderblom et al. 09931 1 for the Pleiades and by Da 
Silva et al. (|2009|) for members of young associations, this may 
just be a result of the selection effects described above. 

To take isuch a dependence on rotation into account, we de- 
rived the mean locus of IC2391 members on the Li EW vs B-V 
diagram and the residuals from it. When plotting such residuals 
versus rotational period, a correlation is found with significant 
results at about 2 cr (Fig.|9]l. HD 61005 fits such relation(s) fairly 
well, i.e. it has the expected lithium content for IC2391 stars with 
its rotational period. A similar trend is also seen considering the 
projected rotational velocity. 

5.3. Cliromosplieric and coronal activity 

There are no members of Argus association with measured H&K 
emission, so no comparison is possible for this indicator. X-ray 
data of IC 2391 are taken from several sources and summarized 
in Appendix A. X-ray data of Argus field members are all taken 
from the ROSAT Bright Source Catalog. Figure [TOl shows the 
distribution of log Lx/Li,„i for members of Argus association, IC 
2391 open cluster, and HD61005. It results that HD61005 has 
a lower X-ray luminosity than most of the Argus and IC 2391 
members, where a few slow rotators have similar coronal lu- 
minosities. A few other members have no X-ray detections in 
ROSAT All Sky Survey. A systematic difference between field 
stars and IC 2391 members is also seen in X-ray emission, most 
likely due to the selection effects previously discussed. 
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5.4. Chemical composition 

Abundance ratios of iron-peak and a-elements for HD 61005 are 
very similar to those derived for IC 2391 members by DR09, re- 
flecting in both cases a solar chemical composition. The mean 
difference in [X/Fe] ratios of our star with IC 2391 is only ~ 
0.02 dex, when considering all the elements listed in Table |4l 
Our analysis was performed in the same way as that presented 
in DR09 and D'Orazi et al. (I2009ab for IC 2391 members, so 
that systematic errors in the abundance scale are not of con- 
cern for such comparison. In Fig. [TT] we show [X/Fe] ratios 
as a function of [Fe/H] for HD 61005 along with the average 
values and their corresponding standard deviation as obtained 
from IC 2391 stars. Similarly, Argus association field stars have 
a nearly solar iron content (Viana Almeida et al. 2009l. This 
provides further support for the link between HD 61005 and the 
Argus association. However, concerning Fe, Na, Si, Ca, Ti, and 
Ni, IC 2391 does not have a peculiar abundance pattern with re- 
spect to nearby young stars: all the young clusters/associations 
surveyed so far from a chemical point of view show a compo- 
sition that is very close-to solar (Santos et al. 120081 D'Orazi et 
al. 12009b ', Biazzo et al. 1201 11 1. This means that the chemistry 
cannot conclusively state the physical association of HD 61005 
to Argus. 

On the other hand, a very interesting datum is given by 
the Ba content: HD 61005 shares with IC 2391 dwarfs the un- 
usual, extremely high [Ba/Fe] ratio, which is about a factor of 
five above the solar value ([Ba/Fe]a!0.7 dex). As discussed in 
D'Orazi et al. (2009a), this anomalous Ba abundance is likely 
due to NLTE effects (probably arising from over-ionization phe- 
nomena caused by the high level of chromopsheric activity that 
is typical of young solar type stars, see e.g. Schuler et al. 1201 Oi l 
rather than to a real higher content of neutron-capture elements 
in young, nearby stars. However, our result points out that the 
[Ba/Fe] trend in HD 61005 is very different from the one derived 
for old, solar-type stars, providing an independent, complemen- 
tary observational hint of the young nature of this object (see 
D'Orazi et al. 20(39a for a detailed discussion on the relationship 
between Ba and cluster age). 

5.5. IR excess 

As mentioned in Sect.[T] HD 61005 has a very large IR excess. 
We briefly discuss this in the context of our revised age of the 
system, because the occurrence of an IR excess and its amount 
with respect to the photosphere are both observed to decline with 
age. IC 2391 was observed at 24jj.m with Spitzer by Siegler et 
al. (120071 . resulting in a frequency of FGK stars with IR excess 
of 31^g^%, which fits the trend toward decreasing IR excess as a 
function of age, as observed for other clusters. No observations 
of IR excess have been published for the field members of the 
Argus association. HD 61005 has a 24 jjin excess {/excess/ fpimt - 
1 .096 + 0.049; Meyer et al. 2008) that is more than for any star in 
IC 2391 included in the Siegler et al. (2007 ) sample, with the F3 
star HD74374=PMM1174just below this value (/excess/ fphot - 
0.98). 

No clear trend toward IR excess with stellar rotation results 
from IC 2391 data: IR excess is found on slow rotators like PMM 
441^ and PMM 4467, but also on the very fast rotator PMM 

PMM 4413 is similar to HD 61005 in several aspects: it has a 
slightly earlier spectral type, slow rotation (P,.„, = 5.05 d), faint X ray 
emission (logLx/L;,,,; = -4.47), Li EW below the mean locus of the 
cluster and presence of 24 fim excess. However, it is a double-lined 
spectroscopic binary with orbital period 90 days (Platais et al. i20071 l. 



1820. Observations of the whole cluster, as well as the field 
members of the association, might provide higher quality statis- 
tics for this kind of analysis. A similar analysis for the Pleiades 
open cluster is also inconclusive: Greaves et al. (2009) note that 
the four solar type stars with 24 fim excess are slow rotators but 
selection effects might be at work because half of fast rotators 
have no excess measurement, because they are projected in a re- 
gion of heavy cirrus (Gorlova et al. |2006). 

6. Conclusion 

We have examined several properties of HD 61005 with the 
goals of a better characterization of the system and a robust 
determination of its age. We determined rotation period, radial 
and projected rotational velocity, chromospheric emission, ef- 
fective temperature, and chemical composition (see Table |6] for 
a summary of target properties). From projected rotational ve- 
locity and rotation period we inferred that the star is seen close 
to edge-on. In that case, there is no indication of any significant 
misalignment with the debris disk. 

Table 6. Summary of parameters of HD 61005 



Parameter 


Value 


Reference 


V 


8.22+0.01 


this paper 


B-V 


0.751+0.003 


Menzies et al. 119901 


V-I 


0.805+0.005 


Menzies et al. l|1990l 


V-Ks 


1.762+0.027 


this paper 


J-K 


0.447+0.036 


2MASS 


n (mas) 


28.29+0.85 


van Leuveen (120071 


My 


5.478+0.068 


this paper 


L/Le 


0.583+0.048 


this paper 


R/Re 


0.840+0.038 


this paper 


Teff 


5500+50 


this paper 


log.? 


4.5+0.2 


this paper 


[Fe/H] 


0.01+0.04 


this paper 


EW Li 6708 (mA) 


171+3 


this paper 


A(Li) 


2.85+0.1 


this paper 


S index 


0.502 


this paper 


log/?M 


-4.310 


this paper 


logLx 


29.3+0.2 


this paper 


logLx/Lto, 


-4.02 


this paper 


V sin i (km/s) 


8.2+0.5 


this paper 


Pro, (d) 


5.04+0.04 


this paper 


Pcyc (y) 


3.4+0.1 


this paper 


^star 


■7-7+13 
' '-15 


this paper 


RV (km/s) 


22.5+0.5 


this paper 


^lalpha (mas/yr) 


-55.71+0.59 


van Leuveen (I2007J 


fideha (mas/yr) 


74.58+0.62 


van Leuveen ( 20071 


U (km/s) 


-22.2+0.6 


this paper 


V(km/s) 


-14.3+0.3 


this paper 


W (km/s) 


-4.1+0.2 


this paper 


age (Myr) 


40 


this paper 



We found that when applying the standard age calibrations 
to the values of several age indicators (lithium, chromospheric 
and coronal emissions, rotation period) HD 61005 results of age 
comparable to that of the Pleiades (120 Myr). However, the kine- 
matic parameters strongly indicate membership in the Argus 
association, which is significantly younger (40 Myr; Torres et 
al. 2008). We then compared the properties of HD 61005 to 
those of Argus association members, including the open cluster 
IC 2391. We found that HD61005 parameters are on the edges 
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but not outside the distribution of IC2391 members. The lithium 
content and coronal emission are similar to those of members of 
IC 2391 of comparable rotational period. HD 61005 also has a 
similar chemical composition to IC 2391. 

HD 61005, therefore, can have the same age as the slowly 
rotating, less active, and (relatively) Li -poor stars in IC 2391. 
The kinematic parameters strongly indicate that this is indeed 
the case. We therefore conclude that the association with Argus 
is very likely and the age of HD 61005 is about 40 Myr. 

The younger age for HD 61005 than currently assumed (90 
Myr, Hines et al. 2007]l has a significant impact on the detec- 
tion limits of planetary companions in direct imaging programs. 
Buenzh et al. (120101) derive limits of about 3 and 5 My at 50 
and 15 AU, respectively, assuming a 90 Myr age. These limits 
become significantly smaller when taking the revised age into 
account. Therefore, only less massive planets may perturb the 
disk and be a likely cause for the peculiar features described in 
Sect.[T] Dedicated modeling of the possible masses and locations 
of planets that might explain the disk features are needed. 

A younger age also helps in explaining the unusually large 
IR excess observed for HD 61005, since debris disk frequency 
and luminosity are known to decline with age (Evans 120091 
Wyatt 2008 ). However, even at an age of 40 Myr, the 24 fim 
excess of HD 61005 still turns out to be larger than that of any 
IC 2391 member for which suitable data are available. The solar 
abundance of HD 61005 is consistent with the lack of significant 
metallicity differences between stars with debris disks and gen- 
eral field stars, at odds with what has been found for the parent 
stars of giant planets (Greaves et al. 2006). 

Finally, from our comparison of HD61005 with the Argus 
association, differences can be noted between field and cluster 
members in Li EW, X-ray luminosity, and possibly in rotational 
period distributions. Argus field stars are on average more active, 
rotate faster, and have more lithium than their IC 2391 cluster 
counterparts. 

A small age difference (5-10 Myr), with field stars being 
slightly younger than cluster members, might explain the lithium 
difference. A vigorous star formation event followed by super- 
nova explosion(s) might have triggered star formation in nearby 
regions, with some delay after the cluster was formed. This 
kind of scenario is envisaged for the Sco-Cen complex as a 
sequence of star formation events, involving Lower Centaurus 
Crux, Upper Centaurus Lupus, Upper Scorpius, TW Hyd asso- 
ciation, e and Cha (Ortega et al. 2009j. 

However, our favorite explanation for these differences is 
a selection effect of Argus field stars all being selected from 
the SACY sample (Torres et al. 120061 ). which only includes 
stars with X-ray counterparts in the ROSAT Blight Star Catalog 
(Voges et al. 1999), while 1C2391 members are selected using 
a variety of membership indicators. Our adopted target list of 
IC 2391 members from Torres et al. (2008 1 includes additional 
selection criteria aimed at better defining the kinematical prop- 
erties of the association (excluding known or suspected spectro- 
scopic binaries unless an orbital solution is available). 

Our results suggest that, at least for the associations be- 
yond 100 pc in the Torres et al. (2008^ list, a selection crite- 
rion based on inclusion in the ROSAT Bright Star Catalog might 
have caused a bias towards the most active and fast rotating 
stars. Therefore, we suggest that other less active, slowly rotat- 
ing members probably still wait to be identified. 

If the different disk lifetime is one of the causes of the dis- 
persion of stellar rotation rates in coeval groups of young stars 
and if planet formation occurs more easily in long-lived disks, 
the selection bias that we identified might have a significant im- 



pact for disk and planet searches. The coiTelation of lithium and 
other age indicators with stellar rotation might also cause the 
age determination for isolated field stars to be biased if such a 
dependency is not taken into account. 
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Table A.l. X-ray emission of IC 2391 members. 



Star Star IRXS V B-V logLx Ref. 

PMM VXR 



7422 




J082845. 5-520523 


10.40 


0.69 


29.95 


1 


7956 




J082952.7-5 14030 


11.50 


0.97 


30.17 


1 


6978 




J083502.7-521339 


12.01 


1.02 


29.52 


2 


2456 




J083543.4-532123 


12.16 


0.92 


29.94 


2 


351 




J083624.5-540101 


10.17 


0.57 


30.26 


1 


4336 


2b 


J083756.4-525714 


11.25 


0.86 


29.98 


2 












29.72 


3 


4362 


3a 




10.91 


0.67 


29.36 


3 


4413 


5 




10.20 


0.67 


29.36 


3 


4467 


12 


J083954.3-525755 


11.80 


0.86 


29.52 


2 












29.63 


3 












29.39 


4 


1083 


14 


J084007.4-533744 


10.38 


0.57 


30.15 


1 












30.17 


3 


8415 


16a 


J0840 16.6-525629 


11.63 


0.87 


30.03 


2 












30.03 


3 












30.04 


4 


1759 


18 




12.54 


1.25 


29.59 


3 


1142 


22a 


J084049.6-533737 


11.04 


0.68 


29.93 


2 












30.26 


3 


1820 


35a 


J084124.0-532247 


12.41 


0.99 


30.22 


1 












30.11 


3 


4636 


41 




13.20 


1.36 


29.66 


3 












29.75 


4 












29.88 


5 


3695 


47 




13.30 


1.44 


29.57 


3 












29.79 


4 












29.72 


5 


2888 


66 




9.76 


0.46 


29.43 


3 


2012 


69a 




11.41 


0.83 


29.85 


3 


4809 


70 


J084407.4-525316 


10.73 


0.63 


29.73 


2 












29.84 


3 


5884 


72 


J084425.2-524219 


11.35 


0.72 


29.96 


2 












30.18 


3 


4902 


76a 


J084526.6-525132 


12.45 


1.04 


29.70 


2 












29.93 


3 


6811 


77a 


J084539.5-522556 


9.91 


0.60 


30.24 


1 












30.18 


3 


X-ray non detections 


1560 






10.64 


0.61 


< 




6974 






12.04 


1.04 


< 




4280 






10.04 


0.67 


< 




3359 






11.38 


0.76 


< 




5376 






13.94 


1.37 


< 




665 






11.33 


0.75 


< 




686 






12.55 


1.04 


< 




756 






11.06 


0.68 


< 




5811 






9.16 


0.37 


< 




1373 






12.06 


0.96 


< 




2182 






10.18 


0.63 


< 





Sources of data: 1: ROSAT Bright Source C atalog (Voges et al. ,1999i l: 
2: ROS AT Faint Source Ca talog ( Voges et al. l2000); 3: Pa tten & Simon 
n996b ; 4: Simon & Patten fT9^ : 5: Marino et al. ([2005}. 



Appendix A: X-ray emission of IC 2391 members 

Because the IC2391 cluster extends over an area larger than the 
one considered in pointed X-ray observations and because sev- 
eral new members are considered, we summarize here the X- 
ray luminosities of the IC2391 members selected from Torres 
et al. (2008 ). We collected the results of pointed observations 
with ROSAT (Patten & Sim on 1996; Simon & Patten 1998) and 
XMM (Marino et al. 120051 ). We also cross-checked the Torres et 
al. (i2008j membership list with the ROSAT all-sky bright and 



faint sources (Voges et al. 119991 120001 ). adopting a matching ra- 
dius of 30 arcsec. X-ray fluxes were derived using the calibration 
by Hiinsch et al. ( 119991 ) and X-ray luminosities adopting the dis- 
tances from Torres et al. (l2008j . 
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Fig. 10. Distribution of log Lx/Ltoi of IC 2391 (filled circles), 
Argus association (blue filled squares), and Pleiades members 
(open squares and open triangles for detections and upper limits 
respectively; data taken from Stauffer et al. 1994) versus B-V 
colour. HD 61005 is plotted as a (red) star. A few IC 2391 mem- 
bers are not detected in X-ray. 



vsini (lim/s) 

Fig. 9. Upper panel: residuals of Li 6708 EW with respect to 
mean locus vs B-V defined by IC 2391 members plotted versus 
rotational period (dashed line in Fig|8]i. Filled and open (black) 
circles: IC 2391 stars with confirmed -i-likely and uncertain rota- 
tion periods, respectively; filled and open (blue) squares: Argus 
field stars, with confirmed -i-likely and uncertain rotation periods, 
respectively; (red) star: HD 61005. The continuous line repre- 
sents the fit of IC2391 members with confirmed rotational peri- 
ods, and the dotted line represents the fit including the uncon- 
firmed periods. Lower panel: residuals of Li 6708 EW with re- 
spect to mean locus vs B-V defined by IC 2391 members plot- 
ted versus projected rotational velocity. Filled circles: IC 2391. 
Filled squares: Argus field. (Red) star: HD 61005. 
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Fig. 11. Abundances of individual elements of HD 61005 com- 
pared to thos e of IC 2391 (from D'Orazi & Randich |2009l and 
D'Orazi et al. l2009al l. 
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13 



Period 


AP 


AV 


Acccur. 


JD range 


#of 


PN @99% 


PN 


Vmin 


Source 


(d) 


(d) 


(mag) 


(mag) 


(2440000.0) 


meas. 


conf. level 




(mag) 




5.04 


0.04 


0.05 


0.033 


11869.8-15137.8 


588 


8.03 


13.1 


8.22 


AS AS full JD range 


5.05 


0.15 


0.03 


0.034 


12684.8-12746.8 


29 


5.71 


6.9 


8.22 


AS AS 


5.08 


0.15 


0.04 


0.034 


13418.7-13479.7 


25 


5.76 


6.0 


8.20 


AS AS 


5.09 


0.15 


0.05 


0.033 


14152.7-14213.7 


22 


5.29 


7.7 


8.19 


AS AS 


5.04 


0.20 


0.05 


0.015 


7857.6-9030.1 


97 


5.56 


10.9 


8.22 


Hipparcos full JD range 


5.02 


0.80 


0.05 


0.014 


8043.8-8254.7 


25 


5.54 


7.1 


8.21 


Hipparcos 


5.07 


0.05 


0.22 


0.108 


7857.6-9030.1 


102 


3.72 


6.6 


8.21 


Tycho full JD range 



Table 2. Measurements of FEROS spectra. OS= Object-i-Sky OC=Object=Calibration (simultaneous thorium). 



JD 


RV 






vsmi 




TelT 






EWLi 


S index 


Remarks 


-2450000 


km/s 






km/s 




K 






mA 






2982.81810 


+22.5 ± 


0.1 


7 


8 + 0.8 


5521 


11 ± 


22 


42 


163. 


0.492 


OS 


2982.82576 


+22.5 ± 


0.1 


7 


8 + 0.8 


5543 


06 ± 


33 


52 


170. 


0.501 


OC 


3447.62607 


+22.5 ± 


0.1 


8 


4 ±0.8 


5496 


60 ± 


22 


20 


172. 


0.518 


OS 


3447.63732 


+22.5 ± 


0.1 


8 


4 + 0.8 


5488 


43 ± 


24 


10 


173. 


0.529 


OC 


3807.55698 


+22.5 + 


0.1 


8 


4 + 0.9 


5487 


95 ± 


26 


98 


173. 


0.520 


OS 


3807.61385 


+22.5 + 


0.1 


8 


5+0.9 


5497 


87 ± 


28 


86 


172. 


0.508 


OC 


3807.69345 


+22.5 + 


0.1 


8 


4 + 0.9 


5485 


21 ± 


28 


06 


171. 


0.506 


OC 


3810.50830 


+22.5 ± 


0.1 


8 


2 + 0.8 


5511 


77 ± 


23 


01 


171. 


0.512 


OC 


3810.58558 


+22.5 ± 


0.1 


8 


2 + 0.9 


5495 


18 ± 


23 


51 


172. 


0.517 


OC 


3810.70741 


+22.5 + 


0.1 


8 


3 ±0.8 


5477 


39 ± 


25 


16 


173. 


0.508 


OC 


3815.56421 


+22.5 + 


0.1 


8 


2 ±0.8 


5481 


98 ± 


24 


40 


173. 


0.507 


OC 


3815.64599 


+22.5 + 


0.1 


8 


2 ±0.8 


5495 


80 ± 


22 


50 


172. 


0.506 


OC 


mean 


22.5 






8.23 




5499 




171.3 


0.510 




rms 


0.0 






0.27 




19 






2.8 


0.010 




night mean 


22.5 






8.21 




5500 




171.1 


0.510 




rms 


0.0 






0.25 




18 






2.6 


0.010 





Table 3. Summary of measurements of spectroscopic parameters of HD 6 1005 available in the literature 



RV 


vsini 


EW Li 


S index 


Instrument 


Resolution 


Reference 








0.492 






Henrvet all 19961 


22.0 


8 


176 




FEROS 


48000 


Wichmannetal.|2003| 


22.3 ± 0.2 


9 






CORAVEL 




Nordstrom et al. 120041 


23.0 


8 






AAT 


51000 


Waite et al.,20051 








0.459 


CTIO 1.5m 


1100 


Gray etal. 120061 


21.98 ±0.85 


14.9 


176 


0.444 ±0.081 


Palomar 


16000 


White etal. 120071 


22.5 ±0.1 








FEROS 


48000 


Setiawanetal. 120081 




8.2 




0.471 


FEROS 


48000 


Schroeder et al.'2009i 




9.9 ±0.9 


169 ±3 




FEROS 


48000 


Weise et al. 2010 


22.5 ±0.1 


8.2 ±0.5 


171 ±3 


0.516 


FEROS 


48000 


this work 


22.5 ±0.1 


8.2 ±0.5 


171 ±3 


0.502 






adopted 



Table 4. Metalhcity and [X/Fe] ratios for HD 61005, where errors represent uncertainties due to EWs. 



Name 


[Fe/H] 


[Na/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Ni/Fe] 


[Ba/Fe] 


HD 61005 
IC 2391" 


0.01±0.04 
-0.01±0.02 


-0.03±0.05 
-0.04±0.05 


0.00±0.04 
0.01±0.02 


0.00±0.04 
0.02±0.01 


0.02±0.06 
0.00±0.02 


0.00±0.06 
0.00±0.02 


0.63±0.06 
0.68±0.07 



" Average abundances (with standard deviations from the mean) for IC 2391 by D'Orazi & Randich [2009l and D'Orazi et al. l2009al for comparison 



